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Cosmological-Model-Parameter Determination from Satellite- Acquired 
Supernova Apparent Magnitude versus Redshift Data 

Silviu Podariu^, Peter Nugent'^, and Bharat Ratra^ 
ABSTRACT 

We examine the constraints that sateUite-acquired supernova apparent magnitude 
versus redshift data will place on cosmological model parameters in models with and 
without a constant or time-variable cosmological constant A. Data which could be 
acquired in the near future will result in tight constraints on these parameters. For 
example, if all other parameters of a spatially-flat model with a constant A are known, 
the supernova data should constrain the non-relativistic matter density parameter 
to better than 1% (2%, 0.5%) at la with neutral (worst case, best case) assumptions 
about data quality. 

Subject headings: cosmology: observation — large-scale structure of the universe — space 
vehicles — supernovae: general 

1. Introduction 

Recent applications of the apparent magnitude versus redshift test based on Type la supernovae 
(SNe la) have resulted in interesting constraints on cosmological- model parameters (see, e.g., Riess 
et al. 1998; Perlmutter et al. 1999; Podariu & Ratra 2000; Waga & Frieman 2000; Gott et al. 
2001). Higher quality data will result in tighter constraints on cosmological- model parameters. A 
dedicated supernova space telescope could provide the high quality data needed to realize the full 
potential of this neoclassical cosmological test. 

In this paper we examine constraints on cosmological-model parameters that will result from 
such a data set. For definiteness we focus on data that could be acquired by the proposed SNAP 
space telescope (Curtis et al. 2000 and http:/ /snap. lbl.gov). That is, we assume a data set of 2000 
SNe la multi-frequency light curves, for SNe out to redshift z = 2, with errors discussed below. 

Observational data favor models with a low ^}q. The simplest such models have either flat 
spatial hypersurfaces and a constant or time- variable cosmological "constant" A (see, e.g., Peebles 
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1984; Peebles & Ratra 1988; Sahni & Starobinsky 2000; Steinhardt 1999; Carroll 2000; Binetruy 
2000), or open spatial hypersurfaces and no A (see, e.g., Gott 1982, 1997; Ratra & Peebles 1994, 
1995; Kamionkowski et al. 1994; Gorski et al. 1998). For a constant A (with density parameter 
Oa), these models lie along the lines + = 1 and = 0, respectively, in the more general 
two-dimensional (Og, model-parameter space. Depending on the values of and models 
in this two-dimensional parameter space have either closed, flat, or open spatial hypersurfaces. In 
this paper we derive constraints on the parameters of the two-dimensional model as well as those 
of the special one-dimensional cases. 

We also derive constraints on the parameters of a spatially-flat model with a time-variable A. 
The only known consistent model for a time-variable A is that based on a scalar fleld (0) with a 
scalar field potential V{(f)) (Ratra & Peebles 1988). In this paper we focus on the favored model 
which at low z has V{4>) oc t^"", a > (Peebles & Ratra 1988; Ratra & Peebles 1988)^. This model 
is in reasonable accord with observational data (see, e.g., Peebles &; Ratra 1988; Ratra & Quillen 
1992; Podariu & Ratra 2000; Waga & Prieman 2000; Brax, Martin, & Riazuelo 2000)^. 

A scalar field is mathematically equivalent to a fluid with a time-dependent speed of sound 
(Ratra 1991), and it may be shown that with V{(f)) oc (f)"'^, a > 0, the ^ energy density behaves 
like a cosmological constant that decreases with time. We emphasize that in our analysis of this 
model here we do not make use of the time-independent equation of state fluid approximation to 
the model that has sometimes been used for such computations (see the discussion in Podariu & 
Ratra 2000). 

Huterer & Turner (1999), Starobinsky (1998), Nakamura & Chiba (1999), Saini ct al. (2000), 
and Chiba & Nakamura (2000) discuss using supernova apparent magnitude versus redshift data 
to determine the scalar fleld potential of the time-variable A model. This is a difficult task. Maor, 
Brustein, &; Steinhardt (2001) note that even data of the quality anticipated from SNAP will 
not result in very tight constraints on an arbitrary equation of state. They consider a simple 
illustrative example, with an equation of state parameter w that has two terms, one constant and 
the other linear in z. Maor et al. show confidence contours (in a two-dimensional plane) for the 
two parameters in the equation of state for this model in their Figure 2. After marginalizing over 
r^o the peak-to-peak spread in their 2 a contour for the equation of state dX z = 0, wq, is about 
—0.3 for Wo = —0.7, or about 43% of the value of wq. This corresponds to a symmetrized 2 a 
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uncertainty of about ±22% on wq^. The corresponding peak-to-peak spread in their 1 a contour 
is about —0.22, which corresponds to a 1 o" uncertainty of about ±16% on wq. For fixed Qq, the 
peak-to-peak spread in their 1 a contour is about —0.09, which corresponds to a 1 o" uncertainty of 
about ±6.5% on wq. While much larger than the constraints we place on model-parameter values 
(see below) this is still a reasonably precise determination of wq. 

Motivated by the approach adopted in analyses of current supernova apparent magnitude 
versus redshift data (see, e.g., Riess et al. 1998; Perlmutter et al. 1999), we instead focus on how 
well future supernova data will constrain parameters of various cosmological models^. 

We want to determine how well supernova data distinguishes between different cosmological- 
model-parameter values. To do this we pick a model and a range of model-parameter values and 
compute the luminosity distance Di,{z) for a grid of model-parameter values that span this range. 
Figure 1 shows examples of Di,{zys computed in the time- variable A model (Peebles & Ratra 1988). 

The error bars on the supernova fluxes are the ones that are most likely to be symmetric (and 

thus allow for the simplest comparison between model predictions and observational data), so we 
work with flux / oc for the comparison between model predictions and anticipated data. 

For our purposes, the constant of proportionality in this relation is unimportant since the SNc in 
the final reduced data set have been made standardized candles (see, e.g., Phillips 1993, and more 
recently Riess et al. 1998; Perlmutter et al. 1999). 

For computational simplicity we assume supernova data from SNAP will be combined to 
provide fluxes and errors on fluxes for 67 uniform bins in redshift, of width Az = 0.03, with the 
first one centered aX z = 0.03 and the last one at z = 2.01. In each bin the statistical and systematic 
errors are combined to give a flux error distribution with standard deviation (7{z). 

To determine how well supernova data will distinguish between different sets of model-parameter 
values, we pick a fiducial set of model-parameter values which give a flux fpiz) and compute 



where the sum runs over the 67 redshift bins and P represents the model parameters, for instance J7o 
and in the general two-dimensional constant A case. Nf^{P) is the number of standard deviations 
the model-parameter set P lies away from that of the fiducial model. This representation (eq. [1]) is 
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exact for the case where the correlated errors between redshift bins for the distance determinations 
are neghgible. 

The error budget is summarized in the next section. Results are presented and discussed in §3 
and we conclude in §4. 

2. Error Budget 

The following provides a brief overview of the constraints that a satellite-based supernova 
program can place on both the statistical and many of the potential systematic errors. For a more 
complete discussion see the SNAP proposal. 

2.1. Statistical Errors 

Currently a single SN la provides a « 16% measurement of the flux 8% in distance) 
(Jha et al. 1999). A large fraction of this uncertainty almost certainly resides in the correction 
for extinction. By going to space one will be able to greatly increase the wavelength coverage 
and precision of the photometric measurements, thereby reducing this uncertainty considerably. 
(Signal-to-noise of ~ 30 could be achieved for a SN at AB(1.0/xm) = 27.0, with systematics in 
the absolute photometry of <1%.) The SNAP satellite has baselined 15 broad-band filters from 
about 0.3 to 1.7 //m in addition to obtaining spectrophotometry near peak for each SN la. A 
conservative estimate of the intrinsic uncertainty for a given SN la with this type of data set would 
be 10% in flux (Ri 5% in distance). There is potential for reducing this even further through 
the identification of additional parameters that constrain the corrected peak luminosity of SNe la 
beyond the single parameter of light-curve shape currently used. Here we will conservatively assume 
that the statistical uncertainty in satellite-based SN la measurements such as these will be 10% in 
flux. y/N statistics on 2000 SNe la over the 67 aforementioned bins would provide an uncertainty 
of < 2% per bin. 

2.2. Systematic Errors 

A major advantage of a space telescope is the much better opportunity for controlling (or 
studying) the many known (and unknown) sources of error. These include environmental effects, 
evolution, intergalactic dust, unusual cases which bias the distribution, etc. See, e.g., Howell, 
Wang, & Wheeler (2000), Aldering, Knop, & Nugent (2000), Croft et al. (2000), Nomoto et al. 
(2000), Barber (2000), Hamuy et al. (2000), Livio (2000), Totani (2000), and Gott et al. (2001) 
for discussions of some of these issues. Without understanding and limiting these sources of error 
an accurate measurement of the cosmological parameters can not be obtained. Here we mention a 
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few of the potential sources of systematic errors and how space-based observations could constrain 
or eliminate them (a more detailed discussion of these and other sources of systematic errors can 
be found in the SNAP proposal). 

Malmquist Bias. This is the sampling bias due to any low-vcrsus-high-redshift difference in 
detection efficiency of intrinsically fainter supernovae. For the aforementioned redshift range, the 
proposed experiment will attempt to detect every supernova in the observed region of sky at 10% 
of its peak brightness, thus eliminating this source of systematic uncertainty. 

Extinction by "Ordinary" Dust. The proposed experiment will attempt to obtain cross- 
wavelength-calibrated data with broad wavelength coverage for each supernova, so that the dimming 
of the spectrum as a function of wavelength can be measured with high signal-to-noise. Further- 
more, SNe la in early-type galaxies with little to no extinction will be targeted to precisely determine 

the intrinsic colors of a SN la at a variety of light-curve shapes (sec Riess et al. 1996 for a study of 
this at low redshift). This would then allow one to study the ratio of selective to total absorption 
from dust and correct for any potential evolution of this ratio as a function of redshift. 

Extinction by "Gray" Dust. It has been suggested by Aguirrc (1999) that certain large (up to 
~0.lA), and possibly needle-like, dust grains can be expelled from galaxies via radiation pressure 
and can have an opacity curve that is shallow in optical bands, thus making them absorptive while 
producing only small color excess. Such dust would lead the unwary cosmologist into underesti- 
mating CIm or overestimating thus producing a systematic bias. If there is gray dust that has 
had insufficient time to diffuse uniformly in intergalactic space, different lines of sight would have 
differing amounts of extinction due to clumping. This would result in an increase of observed su- 
pernova magnitude dispersion, an effect that is not seen in current observations, and could easily be 
detected by a space-based experiment. Furthermore, it is also possible to detect z < 0.5 gray dust 
by comparing optical and near-IR photometry of SNe (both la and II) found in this redshift range 
since the dust is not completely gray and will show a color excess over a large enough wavelength 
range (see, e.g., Riess et al. 2000). 

SN la Evolution. SNe la with different progenitor properties should result in explosions with 
slightly differing properties, even if there is only one mechanism for creating them, and even if 
this mechanism has a set "trigger" such as the Chandrasekhar limit (Hoflich et al. 2000). If these 
differences are not corrected by the light curve width-luminosity relation presently in use, and if the 
distribution of key parameters of the progenitor stars changes with redshift, the SN la explosions 
observed at high redshift could differ in peak luminosity from those at low redshift, leading to a 
systematic error in the determination of the cosmological parameters. However, a dataset acquired 
from a space telescope should allow corrections for these differences, or allow similar SNe la to be 
identified and matched at high and low redshifts thus mitigating against the effects of changing 
progenitor properties. 

One of the wonderful aspects of using SNe la for cosmology is the fact that the supernova 
bares its entire history, from progenitor through explosion, to the observer. Thus the supernova 
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can't hide the effects of evolution since these will make themselves apparent in the light curves and 
spectra. Figure 2 illustrates this statement. This shows the temporal spectral evolution of a typical 
SN la. At very early times one probes the outer, unburned layers left over from the progenitor. As 
seen in Fisher et al. (1997) this epoch displays spectral features from high- velocity carbon left over 
from the original progenitor and could be used to tightly constrain various theoretical models. At 
later times, near peak brightness, we are beginning to probe the layers of the atmosphere which 
show the intermediate-mass elements systhesized in the runaway thermonuclear explosion. Nugent 
ct al. (1995) showed how some of the spectroscopic features of these elements (Si II &: Ca II) nicely 
correlate with the peak brightness of the SN la. Finally, during the nebular phase, we note the 
strong Fe II emission lines at low velocity left over from the radioactive decay of ^^Ni to ^^Co to 
^^Fe. These observations allow one to directly probe the total amount of ^^Ni synthesized during 
the explosion (see Kuchner et al. 1994 and Fisher et al. 1995). 

Curtis et al. (2000) have identified a series of key observable supernova features that reflect 
differences in the underlying physics of the supernova. By measuring all of these features for each 
supernova one should be able to tightly constrain the physical conditions of the explosion, making 
it possible to recognize supcrnovac that have similar initial conditions and/or arise in matching 
galactic environments. The current theoretical models of SN la explosions are not sufficiently 
complete to predict the precise luminosity of each supernova, but they are able to give the rough 
correlations between changes in the physical conditions of the supernovae and the peak luminosity 
(Hofiich, Wheeler, k. Thielemann 1998). These conditions include the velocity of the ejecta (a 
measurement of the kinetic energy of the explosion), the opacity of the inner layers (which affects 
the overall light curve shape), the metallicity of the progenitor (which affects the early spectra), 
^^Ni mass (a measurement of the total luminosity), and ^^Ni distribution (which might lead to small 
effects in the light curve shape at early time). One can therefore give the approximate accuracy 
needed for the measurement of each feature to ensure that the physical condition of each set of 
supernovae is well enough determined so that the range of luminosities for those supernovae is well 
below the systematic uncertainty bound of ~2% when all the constraints are used together (see 
Curtis et al. 2000 for a full description). 



3. Results and Discussion 

For SNAP data cr(z) (eq. [1]) is estimated to be 2% in each redshift bin up to z = 1.7, and 
then increasing linearly with redshift to 10% at z = 2. This is the "neutral" case. The "best" 
case assumes that errors are limited by \/iV statistics (with systematic errors at or below the 1% 
level), giving a{z) = 1% over the whole redshift range. The "worst" case (this is the baseline SNAP 
mission) assumes a{z) = 3% to z = 1.2 and = 10% from z = 1.2 to z = 2. 

Figure 3 illustrates the ability of anticipated space telescope data to constrain cosmological- 
model parameters for the general two-dimensional constant A case. SNAP data with even worst case 
error bars will lead to greatly improved cosmological-parameter determination (see, e.g., Riess et 
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al. 1998, Perlmutter et al. 1999, and Podariu & Ratra 2000 for constraints from current data). We 

note that as expected the contours are elHptical, indicating that one combination of the parameters 
is better constrained than the other orthogonal combination (see, e.g., Goobar &; Perlmutter 1995). 

Figure 4 illustrates the ability of space telescope data to distinguish between a constant and 
a time-variable A in a spatially-flat model. The fiducial model here is a constant A (a = 0) model 
with = 0.28 and 0,\ = 0.72. SNAP data with even worst case error bars will result in greatly 
improved discrimination (see, e.g., Podariu & Ratra 2000 for the current situation). We note again 
that the contours are elliptical. 

Figure 5 illustrates the ability of space telescope data to constrain and a in the spatially- 
flat time- variable A model (Peebles Sz Ratra 1988). Here the time- variable A fiducial model has 
Clo = 0.2 and a = 4. Again, SNAP will allow for tight constraints on these cosmological parameters. 

If other data (such as cosmic microwave background anisotropy measurements from MAP and 
Planck Surveyor and weak-lensing studies from the proposed SNAP mission) pinned down some of 
the cosmological parameters, the supernova data would then be able to provide tighter constraints 
on the remaining parameters. For instance, Figure 6 shows constraints from space telescope data 
on f^o in a spatially-flat constant A model and in an open A = model. As expected from the 
elliptical shape of the contours in Figure 3, anticipated supernova data will constrain Qq more 
tightly in the spatially-flat case than in the open case. In both cases SNAP will provide tight 
constraints on CIq. For instance, at 3 a, in the spatially-flat model we find = 0.3 ± 0.007, 
= 0.3 lb 0.015, and = 0.3 it 0.003 for neutral, worst, and best case errors, while in the open model 
we have fio = 0.3 =b 0.015, = 0.3 ± 0.03, and = 0.3 ± 0.006 for neutral, worst, and best case errors. 

Figure 7 shows the space telescope data constraints on CIq and a in the spatially-flat time- 
variable A model, if other data were to require that either a = 4 or JIq = 0.2. SNAP data will 
provide tight constraints on these parameters. For instance, if a = 4 we find 0,q = 0.2 ib 0.009, 
= 0.2 zb 0.02, and = 0.2 it 0.004 for neutral, worst, and best case errors, while if Oq = 0.2 we have 
a = 4 ± 0.25, = 4 lb 0.5, and = 4 ib 0.1 for neutral, worst, and best case errors, all at 3 a. 

4. Conclusion 

Supernova space telescope data of the quality assumed here will lead to tight constraints on 
cosmological-model parameters. For instance, in a spatially-flat constant A model where all other 

parameters are known, anticipated space telescope supernova data will determine to about 
±0.8%, ±1.7%, and ±0.4% (for neutral, worst, and best case errors respectively) at 1 a. The 
corresponding errors on Oq for the open case are about ±1.6%, ±3.7%, and ±0.7%. For the 
time-variable A model, when a is fixed, will be known to about ±1.5%, ±3.1%, and ±0.7%, 
respectively, while when is fixed, a will be determined to about ±2.2%, ±4.2%, and ±1.2%. 
This will have important consequences for cosmology. 
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Fig. 1. — Lines in the panels in the upper row show luminosity distance Di^{z,a) as a function of 
redshift z for various values of a computed for Hubble parameter Hq = 65kms~^Mpc~^ for the 
spatially-flat time- variable A model with scalar field potential V{(j)) oc In descending order 

at 2; = 2 the lines correspond to a = 0, 2, 4, and 8 (solid, dot-dashed, dashed, and dotted curves 
respectively), a = is the constant A model. Prom left to right the three panels correspond to JIq 
= 0.2, 0.4, and 0.6. The three lower panels show the fractional differences relative to the a = 
case, 1 — Dj^(z, a)/Di^{z, a = 0), as a function of z, for the values of 0,q used in the upper panels. 
Here the lines correspond to a = 8, 4, and 2, in descending order at 2; = 2. 
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Fig. 2. — As the SN la evolves in time the rapid expansion of it's atmosphere allows the observer 
to probe to deeper layers as the optical depth falls off with the diminishing density. At early times 
one views the outermost layers, mostly composed of the unburnt progenitor. Near peak brightness 
the intermediate mass elements of S II and Si II are quite visible. At later times, shortly after 
entering the nebular phase, one views the Fe-peak core of the SN la where the radioactive decay 
of RiO.5 Mq of ^^Ni has taken place since the explosion. 
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Fig. 3. — Contours of = 1, 2, 4, and 8 for the constant A model. Left panel is for anticipated 
supernova data with worst case errors, center panel is for neutral case errors, and right panel is for 
best case errors. The fiducial model is spatially-flat with = 0.28 and JIa = 0.72. 




Fig. 4. — Contours of N^r = 1, 2, 4, and 8 for the spatially-flat time- variable A model (Peebles &; 
Ratra 1988). Left panel is for anticipated supernova data with worst case errors, center panel is 
for neutral case errors, and right panel is for best case errors. The fiducial model has = 0.28 
and Q! = (and is thus a constant A model with Ha = 0.72; this was also the fiducial model used 
for Fig. 2). 
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Fig. 5. — Contours of N^^ = 1, 2, 4, and 8 for the spatially-fiat time- variable A model. Left panel 
is for anticipated supernova data with worst case errors, center panel is for neutral case errors, and 
right panel is for best case errors. The fiducial model has = 0.2 and a = 4. 
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Fig. 6. — N„{0,q) for a flat model with a constant A (left panel) and for an open model with no 
A (right panel). In both cases the fiducial model has = 0.3, with Q,\ = 0.7 and respectively. 
Solid lines are for neutral case SNAP errors while dotted (dashed) lines are for best (worst) case 
ones. 
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Fig. 7. — N„{Qo) (left panel) and N„{a) (right panel) for the spatially- flat time- variable A model 
(Peebles & Ratra 1988). In both cases the fiducial model has Oq = 0.2 and a = 4. Solid lines are 
for neutral case SNAP errors while dotted (dashed) lines are for best (worst) case ones. 



